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Ion transport in biological and synthetic nanochannels is characterized by such phenomena as 
D ■ ion current fluctuations, rectification, and pumping. Recently, it has been shown that the nanofab- 
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transport characteristics 



ricated synthetic pores could be considered as analogous to biological channels with respect to their 

l|, y|. The ion current rectification is analyzed. Ion transport through 

J>^' cylindrical nanopores is described by the Smoluchowski equation. The model is considering the 

r^ ; 
, ^~ <, symmetric nanopore with asymmetric charge distribution. In this model, the current rectification 

in asymmetrically charged nanochannels shows a diode-like shape oi I — V characteristic. It is 

>' 
0^ ' shown that this feature may be induced by the coupling between the degree of asymmetry and the 

o : 

C^ . depth of internal electric potential well. The role of concentration gradient is discussed. 
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I. INTRODUCTION 

The role of ion transport through narrow protein channels for living cells has been widely 
studied. There are many experimental and theoretical attempts in understanding the ac- 
tivity of biological channels at physiological conditions. Among models of cellular electrical 
activity the Goldman-Hodgkin-Katz (GHK) current equation plays an important role. It is 
worth noting that the equation assumes of a constant electric field. There are other mod- 
els in biophysical literature e.g. barrier models with several free-energy barriers within the 
channel J3|, y, |5|. These models are phenomenological, i.e. the potential energy barriers 
used in the models often do not correspond to physical properties of the channel. Although 
they can give good quantitative descriptions of some experimental data p , they fail to fit 
these data sets which include I — V relations measured in asymmetrical solutions [7[ . One 
of the reasons is that their the electrostatics is modeled incorrectly. On the other hand, the 
Poisson-Nernst-Planck (PNP) model, based on the mean field approximation, is an alter- 
native theory. Apart from its own limitations p, Isl, lli| , the 3DPNP model seems to take 
into account electrostatics correctly [7[. However, in general, it is not possible to obtain an 
exact solution to the Nernst-Planck (NP) equation, coupled to the electric field by means of 
the Poisson equation. That is why we consider a simplified system that allows for intuitive 
understanding of the underlying mechanism of examined phenomenon. 

The intention of the paper is to demonstrate that the asymmetry of the internal physical 
field also plays a significant role in permeation. To prove this conjecture we focus on the 
only one source of the asymmetry, viz. the electrostatics. Thus, geometrical effects of the 
channel are purposely omitted by choice of the cylindrical structure of the channel. Chemical 
structure is not considered, either. Factors such as geometry, chemical structure, etc. may 
also be sources of the potential asymmetry. The selectivity filters seem to be governed by 
different mechanisms [ll[, therefore they are not discussed here. We propose a simplified 
continuous model that is able to explain the experimental behavior of interest. 

The studies of the synthetic channels, which in several aspects resemble the biological ones 



12,|l3|, may be helpful in investigating the mechanism of ionic transport. The ionic currents 
through these nanochannels exhibit several peculiarities. One of them is the asymmetric 
and non-linear shape of current-voltage characteristic both for different biological channels 
QQQ and fo. asy.,ne.nca., sKaped s,v,.Ke«c ones HQHQ. TKese a.yn,n,et.e. 



in the I-V characteristics are related, among others, to the nanopumping mechanism 



and to the asymmetry of nanodiffusion 



Q 



2l|. It is worth mentioning that the role of 



spatial/time asymmetry in large number of physical processes was discussed in J2j], and that 
the experiment has shown that rectification occurs for colloidal particles in a microfabricated 
channel with a topological ratchet-like polarity J23|. The above-mentioned results point out 
that the theoretical explanation of current rectification in nanochannels should enclose the 
intrinsic asymmetry of channels. 

The aim of this work is to stress the role of the asymmetry of potential well inside the 
channel in the rectification process and in the diode-like shape of / — V^ dependence. In 
our model the total current / through the channel is driven both by an external and an 
internal electrical field. We shall investigate the shape oi I — V curve for symmetric and 
asymmetric cases of the potential well and for different depths of the well. For this purpose 
we need the given shape of the potential well inside the channel. Thus, the first problem 
we shall consider is that of determining the electric potential (pint inside the channel due 
to prescribed surface density a. For simplicity, we shall consider a coaxial infinitely long 
cylindrical channel and we take the radius R of the cylinder small enough to allow the ions 
go only through the channel along z-axis, which implies the lack of electrolytic solution 
inside the channel (no screening). If we assume the equilibration in the transverse direction 



of the channel f24^ we can operate with the ^-dependent potential (plnt'i^) '^^ ^^^ z-axis 
only. Therefore our problem is reduced from the three-dimensional to the one-dimensional 
description. We shall discuss the evolution of the probability density of the finding an ion 
inside the channel, therefore we shall use the 1-D Smoluchowski equation. 

The Smoluchowski equation gives the conditional probability that the particle starting 



from the point z{to) reaches the point z at the time t [25|,|26|. It describes the diffusion of 
probability, because the process of diffusion is the superposition of Brownian motions of the 
molecules of the substance under consideration [2a, |22| . In an open state, the measurable 
electric current through the biological channels is at the picoampere level which gives about 
10^ ions per second |2S|- Thus one ion passes the channel at the tens of nanosecond and after 
10^ passages we get the measurable electric current. These estimations enable us to pass 
from one-ion description via probability density, i.e. from the Smoluchowski equation to the 
continuous description in terms of the electric current density, i.e. to the Smoluchowski- 
Nernst-Planck equation. 



II. MODEL 

The problem which we want to consider is that of electric current flowing through cylin- 
drical channel. There is a broad collection of papers in which the cylindrical synthetic 



nanotubes 



mm, 



3C 



3l| and biological cylindrical ionic channels, e.g. gramicidin channel 



Class 1 porins J33| are discussed. Let us consider a dielectric membrane (biological 
or synthetic) separating two large regions of space that contain some electrolytic solutions. 
There is a pore in the membrane (one-dimensional channel of a radius R and a length L, 
see Fig. 1) through which the ions can move more or less freely. 




FIG. 1: Channel in cartesian coordinates {x,y,z) with the z axis along the axis of the cylinder; 
the cylinder has a radius R and a length L; Vq, Vl are potentials and cq,cl concentrations on the 
left and on the right side of the channel, respectively. 



The biological ion channels are ion selective [3J|. The same feature is displayed in some 
synthetic channels |2,|35|- Therefore, for simplicity, we assume that only positive ions can 
enter our pore. In Cole's model 36[ of electrical activity in membranes the driving force for 
the ions is given by the difference between the membrane potential U = \^„(0) — Vout{L) 
and the reversal potential, which is the potential at which the current is zero. In the 
Goldman-Hodgkin-Katz model (GHK model) of ion flux through the channel the sinaplifying 
approximation is made that the potential gradient through the channel is constant [37[ . On 
the other hand in experiments with synthetic membranes the driving force is caused by 
electrodes of potential V^ and V^_, located in the electrolyte at macroscopic distances from 
the membrane (Fig. 1 in |l8|). Thus in both cases we put the related potential in the form: 



W(^) = V+~ Uz/L, 



(1) 



where t/ = V+ — VC_ and L is a length of the channel. (In general, ions moving through the 
channel affect the local electric field. Thus, the related electric potential may be given by a 
more complicated function. However, if the channel is short or the ionic concentrations on 
either side of the membrane are small this approximation seems to be correct p.) 

We consider charged channel with a given shape of internal potential (plnt'iz) along the 
z axis. Therefore, the total electric field E acting on ions inside the channel is a linear 
superposition of the external Eext and the internal Eint fields, the latter one induced by 
charges on the wall of the channel. The goal of the following analysis is to demonstrate that 
the field asjTiimetry is sufficient to observe diode- like shape of J — ^ relation. Let us consider 
two shapes of the internal potential (p^nt^{z) described below and presented in Fig. 2. They 
roughly correspond to two types of channels, short biological (Fig. 2a) and long synthetic 
ones (Fig. 2b). 
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FIG. 2: The shape of potential well (p^^^ (we use diinensionless coordinate z = z/L and potential 
0j.„t = Pe(l)int a) for Cout = 2eo, eint = 80eo, dg = 0,015 e/nm^, R = 0.004 and g = 0.5 ((l)-line) 
and 'g = 0.25 ((2)-line); b) for ratchet-like shape with (f){Li = 0) = (j){L2 = 1) = 0, (j){g) = —1 and 
g = 0.5 ((l)-line) and g = 0.008 ((2)-line). 



The first potential shape we propose to consider is the potential (pint of the form shown 
in Fig. 2a. This shape of the potential well, which results from charged residues localized at 
z = g, might simulate the situation in short biological [e.g. 33[] and synthetic J3l| channels. 
To find this shape we use the Eq. IA.5I (see Appendix) and we put the surface charge density 
into the formula ()A.7|) to be a{z) = (yg6{z — g). The value of the parameter g gives us the 
asymmetry of the function (pint- 



The second shape, which corresponds to continuous charge distribution, is the electro- 
static potential (pint in the "ratchet-like" shape (Fig. 2b). This case seems to be more 
adequate to the long synthetic tubes in which the ratio of 
the approximation of an infinitely long cylinder works well 



ength to width goes to 0, thus 
3. 



&{z) 



9 — ^1 

0(j) A^"(. - L.) Ar = '"^"'f^'''^ for z € (g, U] (2) 

elsewhere 



where g is the asymmetry parameter of this function, (j){g) is the value of the potential for 
z = g. The asymmetry of the potential is controlled by the value of the parameter g (for 
g = (Li + L2)/2 and Li = L — L2 we have the symmetric function, the most asymmetric 
case we have ii g -^ Li (or g —^ L2)). 

III. RESULTS 

The kinetic Smoluchowski equation is commonly used in various physical, chemical, etc.. 



problems [25[. The equation describes diffusion in any physical field. Therefore a model 
based on this equation can be placed between these two GHK and PNP models (the GHK 
equation is the ID Smoluchowski equation in a constant field). 

Let us start from the ID Smoluchowski equation that contain the electrostatic field: 

d ^' 

j(z,t) = -D—p{z,t)+nE,p{z,t) 

d 
where E^ = — — 0(z) is the 2;-component of electric field, j{z,t) denotes probability density 

current, which describes the fiux produced by diffusion of cations. They are driven by 

the difference in their probabilities of entering the pore from the left p{0) and from the 

right p{L), respectively and by ionic migration caused by difference in electric potential. 

We assume here that the mobility of a particle u fulfills the Nernst-Einstein equation, the 

diffusion coefficient D = u/ (3, where (3 = l/ksT. The ion mobility /i = eZcU, where e is the 

electric charge, and Zc is the cation valence. 



We can write the Eq. El in more convenient form: 

MlA . |fl,-./«(..|e-«"V(.,t) = -^j(z,t). (4) 

The measurable quantity is the stationary mass current Ji {Ji ^ Ji{t)) flowing through 
the channel cross-section A = vri?^: 

J, = TcR^jiz) (5) 

(we use here the identification mentioned in the introduction i.e. one ion passes the channel 
at the tens of nanosecond and after 10® passages we get the measurable electric current). 

For simplicity we use dimensionless coordinate 'z = z/L, and potential (f) = /3e0. For the 
potassium cations Zc = 1- Thus, the electric current resulting from mass current is: 

I = FJ = FnR'j-^—^ 1 , (6) 

Jo 

where F is the Faraday's constant. If we put (p^^^ = we obtain the famous Goldman- 
Hodgkin-Katz (GHK) current equation. 

In the case of ratchet-like shape of potential 0j„j Eq. |21 with Li = and L2 = 1 we get 
the analytic solution: 



I = FJ = 

n -, , -, , - /I _p0{5)-^{O)-Fg p^(9)-^(l)~C7g_ -U 



L V u-m^<p u-m^<p' 



(7) 



Now, we analyze I — V relations for two kinds of potential function 4>int{z) that are 
described in Sec. Model. Let us start from the second case being ratchet-like function of 
z. We use the analytic solution for current that is given by Eq. [7| where we put 0j„j(Li = 
0) = 0i„i(L2 = 1) = 0, ;R = R/L = 0.0004, and Dk+ =2* lO^^ mVs. The ion current 
through the nanochannel with no gradient of concentration (cq = Ci = 0.01 M) depends on 
both the depth of the potential well and on the degree of asymmetry of internal potential 
(pint (see Figures 3). For the fixed depth of potential well 0(^) = —9 and different values of 
'g (from 0.5 to 0.00001) we get various degrees of rectification. It can be clearly seen that 
with increasing asymmetry {'g — > 0) the I ~ V characteristic tends to a diode-like shape 
(Fig. 3a). However the depth of potential well (f){g) is important as well. For the lower value 



of 0(^) = —1 (Fig. 3b) we observe non-linear shape oi I — V curve for both g = 0.5 and 
'g = 0.00001. For the asymmetric potential well the I — V characteristic shows asymmetry. 
On the other hand, for the higher value of (j){'g) = —9 (Fig. 3c) the difference between 
symmetric and asymmetric case is much stronger. In Fig. 3c the I — V curve clearly shows 
the rectification (diode-like shape). The direction of the rectification depends on the value 
of g i.e. if it changes from 0.5 to 0.00001 or from 0.5 to 0.99999. 
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FIG. 3: / — y dependence for the (/)j„j being ratchet-like function of z, for 0j„j(Li = 0) = (^^^^{L2 = 
1) = 0, i? = R/L = 0.0004, CO = ci = 0.01 M, Dk+ =2* 10"^ mVs; a) lfi{g) = -9 and g = 0.5 
- solid line, 'g = 0.3 - dashed line, g = 0.15 - dotted line, 'g = 0.00001 - dot-dashed line; and 
b) (pig) = —1, c) (pig) = —9; in both cases for g = 0.5 - symmetric shape of potential function 
((l)-line) and for g = 0.00001 - the asymmetric shape ((2)-line). 



For the potential 0j„( showed in Fig. 2a (where R = R/L = 0.004), the values of 
boundary condition 0(0) and 0(1) differ from 0. What is more in the asymmetric case ((2)- 
line in Fig. 4) the difference between them grows with increasing value of the potential well 
depth 0(^) i.e. for 0(^) = -10, 0(0) = -0.608 and 0(1) = -0.203 whereas for 0(^) = -40, 
0(0) = —2.430 and 0(1) = —0.811. For cq = ci and symmetric potential function 0j„( the 
I — V dependence is symmetric and weakly non-linear in both cases of 0(^) ((l)-line in Fig. 
4). For the asymmetric potential 0j„(, we observe a weak non-linearity and the appearance 
of the reversal potential for the channel (Fig. 4a) which for increasing value of (j){g) grows 



up (Fig. 4b). Therefore in that case we do not obtain the diode-hke shape for the I — V 
curve. Note that these results reproduce the data reported in 
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FIG. 4: I — V curve for the (pj^^^ being localized function of z, for R = R/L = 0.004, cq = ci = 
c = 0.01 M, -Dx+ = 2 * 10~^ m^/s, the depth of potential well and the boundary conditions: a) 
0(^) = -10, 0(0) = ^(1) = -0.302 for (l)-line, 0(0) = -0.608, 0(1) = -0.203 for (2)-line; b) 
0(^) = -40, 0(0) = 0(1) = -1.208 for (l)-line, 0(0) = -2.430, 0(1) = -0.811 for (2)-line; in 
both cases for g = 0.5 - symmetric shape of potential function ((l)-line) and for g = 0.25 - the 
asymmetric shape ((2)-line). 

In described-above cases we put cq = Ci. However the boundary condition 0(0) 7^ 0(1) 
simulates the gradient of concentration cq 7^ Ci in Eq. IHl Thus putting Ci/cq = e'^^^y e'^^^'' 
recovers the diode-like shape (Fig. 5). 
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FIG. 5: I — V curve for the 



'-'int 



being localized function ol z, g = 0.25 - the asymmetric shape 



of potential function, R = R/L = 0.004, cq = 0.01 M, D^+ = 2 * lO"-* mVs, (l)-line for ci 



Co, 



(2)-line for ci = cqc'^^^' / e'^^^' , the depth of potential well and boundary condition: a) 0((7) = —10, 
0(0) = -0.608, 0(1) = -0.203; b) 0(5) = -40, 0(0) = -2.430, 0(1) = -0.811. 



One can obtain similar results of / — V^ dependence for several other shapes of (pmt with 
cases of boundary conditions discussed-above. 
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FIG. 6: / — y dependence for the (/)j„j being ratchet-like function of z, for <^^^^{L\ = 0) = (/)j„4(L2 = 
1) = 0, i! = RjL = 0.0004, Co = ci = 0.01 M, 0^+ = 2* lO^^ mVs and g = 0.00001; a) lp{g) = -5; 
b) 0(5) = -50; c) 0(5) = 5; d) ^{g) = 50. 



Finally, let us analyze the difference between the potential well and the potential barrier 
inside the channel. The cation that passes the negatively charged channel experiences po- 
tential well, in the contrast the anion goes over the potential barrier. We can see that the 
direction of the rectification is different in these two cases (Fig. 6). What is more important, 
the current increases with the growth of the depth of the potential well, whereas the growth 
of the height of the potential barrier results in dramatic decrease of the current (Fig. 6). 
The main point is that these two cases (potential well and barrier) are not fully symmetric. 
The contrast between these two cases leads us to the conclusion that the same electric field 
could cause enhancement of the cation current and inhibition of the anion current therefore 
the resulting cation selectivity of the channel. However, it has been shown recently that 
the selectivity of ionic channels may originate from physical mechanisms different from the 
electrostatics interactions, too (ll|. 
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IV. COMMENTS AND CONCLUSIONS 

We found that the asymmetry of the internal potential 0j„t coupled with the potential well 
depth are sufficient for the rectification observed as a diode-like shape oi I — V characteristic. 
The concentration gradient is not required for that effect. However, in the absence of internal 
potential 4>int{z)), the Goldman- Hodgkin-Katz current equation, Eq. IHl gives ohmic curve 
I — V a Ac = and non-ohmic shape if Ac 7^ J37|. In the latter case the deviation from 
ohmic behavior results from the gradient of concentration on both sides of the channel and 
is weakly non-linear. In our model, when we put Ac = 0, the internal potential 0j„j (with 
the boundary values 0j„i(O) = 0j„t(l) = 0) is the source of asymmetry, and gives strong 
non-linearity of / — V^. Note that by varying these two factors the degree of the rectification 
can be intensified or diminished. 

The conclusions are supported by the fact that we solve the Laplace's equation instead of 
the nonlinear Poisson-Boltzmann equation (see Appendix). It allows us to avoid a misleading 
intuition that a nonlinear diode-like shape oi I — V characteristics is caused by the non- 
linearity of Poisson-Boltzmann equation. Apart from that, in general, it is not possible to 
obtain an exact solution to the Nernst-Planck (NP) equation, coupled to the electric field 
by means of the Poisson-Boltzmann equation. The purpose of this paper is to show that 
a simplified continuous model that makes use only of the electrical properties of the open 
channel may support our intuitive knowledge of examined phenomenon. 

Moreover, our understanding of the mechanisms which account for experimental obser- 
vations will enable us to design both artificial biological and synthetic channels with desired 
properties. Engineered nanopores may have significant applications. The properties of 
channels and pores that might be engineered include conductance, ion selectivity, gating 
and rectification, and inhibition by blockers 38[- This prospective additionally confirms a 
need for the theoretical models of nanopores. 
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APPENDIX 

The cylinder has a radius R and a height L, the top and bottom surfaces being at z = L 
and z = 0. The surface charge density on the side of the cyhnder is a and we assume that is 
symmetric about the axis of the cyhnder, a = cr(z). We want to find potential at any point 
inside the cylinder. 




FIG. 7: Channel in cartesian coordinates (x, y, z) with the z axis lying along the axis of the cylinder; 
the cylinder has a radius R and a length L; eout-, ^mt are dielectric constants outside and inside the 
channel, respectively. 



There is no loss of generality if we take the infinitely long cylinder charged to prescribed 
charge density c, which can differ from only for z G [0, L] (Fig. 7). In cylindrical coordi- 
nates (r, ^, z) the Laplace's equation takes form 



r 92 Id 1 92 d 

Q^2 J, Qj, ^2 QQ2 Q£2 



0. 



The separation of variables is accomplished by the substitution 
leads to the two ordinary differential equations: 

d'^P IdP , 2 „ 



[r,z) 



(A.l) 
P{r)Z{z). This 



dZ 

dz^ 



(A.2) 



+ k^Z = 0. 



We denote by (pint the potential inside the channel and by (pout the potential outside the 
channel: 



(pout{r, z) = /+^ dk[C{k)e^>^^ + D{k)e-^'^^]K^{\k\r). 



(A.3) 
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Further, these functions must satisfy the boundary conditions 

(f^intir, z) = (f)out{r, z) on r = R 

eo«t ^ emt ^ = -o-(^) on r = R, 

where eouu ^int are the dielectric constants outside and inside the channel, respectively. 
The result on the z axis is 

0Lt(^) = -^^ / dke^'^^a{k)W{\k\), (A.5) 



where 



and 



'^^1 1^ \k\[e,r.M\k\R)Koi\k\R) + e,^M\k\R)Kii\k\R)] ^ '^^ 

1 /'+00 

a{k) = -^= / rf^a(2)e-*'=^ (A.7) 

v27r 7_oo 
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